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Whereas boson coherent states with complex parametrization provide an elegant, and intuitive 
representation, there is no counterpart for fermions using complex parametrization. However, a 
complex parametrization provides a valuable way to describe amplitude and phase of a coherent 
. . . beam. Thus we pose the question of whether a fermionic beam can be described, even approximately, 

' by a complex-parametrized coherent state and define, in a natural way, approximate complex- 

' parametrized fermion coherent states. Then we identify four appealing properties of boson coherent 

, states (eigenstate of annihilation operator, displaced vacuum state, preservation of product states 

under linear coupling, and factorization of correlators) and show that these approximate complex 
(-H ' fermion coherent states fail all four criteria. The inapplicability of complex parametrization supports 

, the use of Grassman algebras as an appropriate alternative. 
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^ ■ I. INTRODUCTION 

> 



As fermion optics research grows, spurred for example by applications to fermionic quantum information process- 
ing the formal question of how to deal with fermionic coherence mathematically becomes important. For bosons. 



J — ' the natural way to describe coherence and optics is in terms of complex-parametrized coherent states, which are 
. eigenstates of the mode annihilation operator, and then employ the optical equivalence theorem A similar 

' analysis for fermions is forbidden because the only complex-parametrized eigenstate of fermion annihilation opera- 
I tors is the vacuum state, but formally a well-behaved fermion coherent state can be introduced by parametrizing 
' with Grassman numbers [6j rather than complex numbers, which overcomes challenges due to anticommutativity rela- 
O il tions 0- These 'Grassman coherent states' are eigenstates of the fermion annihilation operators with Grassman- valued 
eigenvalues and formally demonstrate many of the desirable coherence properties of bosons Q . 
' Although the Grassmann coherent state representation provides a beautiful analogy to the boson coherent state and 
^ , its use in quantum optics experiments , the disadvantage is that the Grassmann algebra and the resulting Grassmann 
coherent state do not convey a clear physical picture of coherence. Complex parameters can be understood in terms of 
amplitude (the modulus of the complex parameter) and phase (the argument of the complex parameter) of the field. 
Here we explore whether complex parametrization can be employed in some advantageous way to analyze coherence 
of fermionic systems. This analysis entails defining an approximate complex-parametrized fermionic coherent state, 
identifying the four chief desirable properties of coherent states, and show that these approximate fermionic coherent 
' states dismally fail every single criterion, which leads to the conclusion that complex parametrization is not useful 
even for approximations to fermion coherent states for any of the four criteria. This strong and negative result further 
reinforces the importance of using Grassman numbers to connect fermion and boson coherence in one formalism. 

Our paper is organized as follows. We review the number state representation for bosons and fermions in Sec. ITTl 
and review boson coherent states in Sec. IIIII Then, in Sec. IIVI we show why complex fermion coherent states in 
a direct analogy with the boson case are impossible and we make an attempt to introduce "approximate fermion 
coherent states" . In Sec. [V] we derive a number of properties of normally-ordered fermion correlators and based on 
these properties we show that no multi-particle fermion state can, even approximately, factorize correlation functions. 
Finally we present our conclusions in Sec. IVII 
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II. NUMBER STATE REPRESENTATION 



111 this paper we denote the boson and fcrmion annihilation and creation operators by a, and c,&, respectively. 
The boson (fermion) operators satisfy the following commutation (anticommutation) relations 

[ai,a]] = Sijt, [ai,aj] 0, [a|,a]] = (1) 
{£,;,£]} = {c,;,c,} = 0, {cl,c]} = (2) 

with 5ij the Kronecker delta and i,j labeling an orthonormal set of modes. 

Number states are states with a definite number of particles in a given mode. They are eigenstates of the particle 
number operator [n = d^a for bosons and n = c^c for fermions) such that n\n) ~ n\n) . The lowest number state, 
the vacuum |0), has no particle, and higher number states can be generated from |0) by the action of the creation 
operator that raises the particle number by one. For bosons, this action is as follows, 

dV) = VnH~T|n+ 1), n = 0,l,2,... (3) 

while for fermions 

ct|0) = |l), ct|l) = (ct)2|0)=0 (4) 

and the latter formula follows from the anticommutation relations 10). Hence, for bosons the number of particles n 
can have any non-negative integer value while for fermions the only allowed values are and 1. This reflects Pauli's 
exclusion principle, which states that there cannot be two fermions in the same state. The action of the annihilation 
operator on number states is the following: 

a|0) = 0, a\n) = y/^\n-\) (n=l,2,...) 

c|0)=0, c|l) = |0). (5) 

The extension of single-mode number to multi-mode number states in general utilizes an arbitrary, but established, 
ordered set of single-mode states |{n}) with ordered occupation numbers {n} = {ni, n2, ns, • • • }. The established 
ordering accommodates the potential sign changes that occur due to the particular algebra chosen to represent the 
quantum statistics of the system. In practice, this ordering is determined through an arbitrary, but established, 
ordering of the creation and annihilation operators which act on the vacuum state |0) (no ordering required now 
within the kct), and the consequence of this ordering is maintained through the use of the appropriate operator 
algebra used when manipulating the creation and annihilation operators. 

In the case of multi-mode boson fields, the commuting algebra introduces a simplification: the operator ordering 
for different modes becomes irrelevant for commuting operators. For example, the state with one boson in mode 1 
and one boson in mode 2 can be obtained from the two-mode vacuum in two equivalent ways: 

|lil2> -dla^lOiOa) (6) 

|l2li> =4al|0i02) (7) 
|lil2> = II2I1). (8) 

Although one may consider the ordering in Eq. ^ to be the established ordering, there is no distinction between the 
resultant state in Eq. © and the one in Eq. lO which uses a different ordering. This is a direct consequence of the 
commuting algebra used for boson fields. It is for this reason that one may use the number state notation [left-hand 
sides of Eqs. © and Q] and the operator notation [right-hand sides of Eqs. © and ^] interchangeably with boson 
fields without regard for the established ordering between different modes. 

However, ordering is important for fields which do not use a commuting algebra. For example, the state with one 
fermion in mode 1 and one fcrmion in mode 2 can be obtained from the two-mode vacuum in two different ways: 

II1I2) =cI4|0i02) (9) 
II2I1) =4cl|0i02) (10) 
II1I2) = -II2I1). (11) 

In the fermion case, the resulting states are not identical, II1I2) 7^ II2I1), because of the fermion anti-commutation 
algebra {cl,^} = (see Eq. (|2l). This illustrates the point that in general an arbitrary, but established, a priori 
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ordering is required. From the practical point of view, it is better to identify a multi-mode number state by a product 
of creation operators that operate on an unordered vacuum state (see the right-hand sides of Eqs. ©, UlUfl ') rather 
than by the ordered occupation- number notation (see the left-hand sides of Eqs. Q and COl)- 

The question of ordering fermion modes explained above becomes unimportant for mixtures of fermion number 
states, as we show in Appendix^ This is significant as it is mixtures of number states, usually in the form of chaotic 
states, that are most often encountered in a real physical system. Chaotic states of bosons and fermions are discussed 
in more detail in Appendix IdI 



III. BOSON COHERENT STATES 



Boson coherent states can be introduced is several ways that emphasize its different physical or mathematical 
properties. Here we consider four interrelated ways. 

Eigenstates of the annihilation operator: — A common definition is that the boson coherent state is an eigenstate 
of the annihilation operator, 

a\a) = a\a), a G C. (12) 
This yields the expansion of a coherent state in terms of number states: 

CXD „ 

To define a multi-mode boson coherent state, consider the complete set of boson modes indexed by k that can be 
assumed to form a countable set fc G N. A multi-mode boson coherent state is then an eigenstate of all annihilation 
operators a^. 

Displaced vacuum state: — The boson coherent state can also be defined as a displaced vacuum state 

|a) ==£)(a)|0) =e""*-"*'^|0). (14) 

In other words, the boson coherent state is a member of the orbit of the vacuum state under the action of the 
Heisenberg-Weyl group D{a), which is parametrized over the complex field. This definition readily extends to multi- 
mode fields where Dk{(y) is the group action for the fc*'^ mode, and the displacement operators for different modes 
commute with each other. 

Linear mode coupler: — The third option for defining boson coherent states is connected to their behavior under 
the action of a linear mode coupler. A linear mode coupler is an important element for transforming both boson 
and fermion fields and it is realized by a beam splitter for optical fields and as a quantum point contact for confined 
electron gases. Furthermore linear loss mechanisms can be described as linear coupling between the signal mode and 
the environmental modes. For bosons, coherent states remain coherent states under linear loss. 

The boson coherent states have the following special property. They are the only states, when mixed on a linear 
mode coupler with the vacuum state, yield a two-mode product state, which is at the same time the two-mode product 
coherent state [T^ . This property is important because it tells that multi-mode coherent states remain as multi-mode 
coherent states under the action of mode coupling and linear loss. 

Factorization of correlators: — The last possibility considered here for defining boson coherent states is related to 
the normally-ordered correlators (correlation functions) 01 

G^''\xi,...,Xn,yn,---,yi) = (V'^(a;i)---V)t(a;„)Vi(y„)-->(2/i)), (15) 

whose complete factorization is identified with the property of coherence. Here ipi^) is the boson annihilation operator 
at the point x, i.e., the sum of annihilation operators of the complete set of modes weighted by the spatial mode 
functions. 

A boson coherent state is then a state for which all the normally-ordered correlators factorize, i.e., for which 

G^"\xi, . . .,xn,yn, . . . ,yi) = f{xi) ■ ■■ f{xn)g{yi) ■ • •g(y„). (16) 

This definition is especially important as the concept of boson coherent states was developed by Glauber 0, ^| to 
describe coherence of the electromagnetic field, and coherent states were defined to factorize the normally-ordered 
correlators in analogy to classical states of the filed. 

In fact there is some ambiguity in the definition of the boson coherent state according to this factorizability condition. 
With the displaced vacuum state in Eq. (|14|) . the phase of the coefficient for Fock state \n) is narga, but the correlator 
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factorizes for an arbitrary phase while the magnitudes of the coefficients have to remain the same. Generahzed 
coherent states that satisfy the coherence condition were studied by Titulaer and Glauber The boson coherent 
state (|14|l is a special case of states in this class. 

Summary: — Each of the four definitions given above are equivalent, which we prove in Appendix In the 
following section we consider extending each of these approaches to constructing coherent states to the fermion case. 

For completeness, we mention another appealing feature of boson coherent states, namely the Gaussian statistics 
governing field quadrature measurements, for example by homodync detection |l2l Il3l|. For fermions, there is no 
obvious way to extend the notion of quadratures so we omit this approach to constructing fermion coherent states. 



IV. FERMION ANALOGY TO THE BOSON COHERENT STATE 



Eigenstates of the annihilation operator: — For fermions, the only eigenstate of the annihilation operator with 
complex coefficient is the vacuum state |0). Indeed, acting on a general pure state = 7|0) + (5|1) by the annihilation 
operator yields the resultant state 5\Q), which is not a multiple of unless 6 = 0. However, if |i5| is small, then c\ip) 
is close to some multiple of \ip), which motivates the consideration of "approximate complex fermion coherent states" . 

One may define an e-complex fermion coherent state as follows. For a given < e < f, an e-complex fermion 
coherent state is a normalized state \u)^, which satisfies 

min c|M)e — < e. (17) 



Here 1 1 IV") 1 1 denotes the norm of a vector, namely 1 1 IV") 1 1 = \'vfW)- 
Proposition 1. For e < 1, the normalized state 

\a) = (cos |a| + sin \a\ e'''''e"c^)|0) (18) 
is an e-complex fermion coherent state i/sin^ \a\ < e. 

Proof. Minimizing the norm of the vector c\a) — f3\a) over /3 € C, we find that the minimum occurs for (3 ^ f3Q ^ 
i sin2|Q;|e'^''S". Evaluating then the norm for this particular Pq, we obtain ||c|q;) — /3o|q!)|| = sin^laj. Hence, if 
sin^ \a\ < e, then \a) is an e-fermion coherent state. □ 

If e 1, one may refer to the e-complex fermion coherent state as an "approximate complex fermion coherent 
state" ^3- Clearly, e bounds the occupation probability for a single particle, and in the limit e — > the set of 
"approximate complex fermion coherent states" reduces to the vacuum state. 

However, attempting to generalize such a definition of the "approximate complex fermion coherent state" to a 
many-mode system introduces important difficulties. Because the fermion annihilation operators of different modes 
are non-commuting, it is in principle not possible to find a common eigenstate, with complex eigenvalue, even if 
a non-trivial eigenstate of the annihilation operator of a single mode existed. Further, to generalize Eq. H18|l . one 
would need to prescribe the ordering of the actions of the operators on the vacuum, which is a step without physical 
justification. This problem will be discussed in the next paragraph. 

Displaced vacuum state: — It is easy to show that Eq. 1)18(1 can alternatively be expressed as 

|a) =e"^'-"*^|0),aGC. (19) 

We note that exp(ac^ — a*c) is formally identical to the boson displacement operator [see Eq. (|14|l ]. We can therefore 
call D{a) = exp(Q!ct — a*c) the complex fermion displacement operator. Any single-mode complex fermion pure state 
can be obtained from the vacuum state by applying a particular D{a): that is, every pure fermion state is in the orbit 
of the vacuum under the action of the displacement operator. If fermion coherent states were defined as displaced 
vacuum states, then any single-mode fermion pure state would be a coherent state, which would not yield a very 
useful definition of coherent states. However, one may still define "approximate complex fermion coherent states" 
according to Eq. ((T^ with \a\ <^ 1. 

To generalize this definition to a multi-mode system, one would have to be careful. A naive approach would be to act 
consecutively by displacement operators of the individual modes on the vacuum, without considering the ordering of 
these operators. However, displacement operators for different modes neither commute nor anticommute so ordering 
infiuences the resultant state. For a given set of ai, q;2, ■ ■ ■ , and a given permutation P of the modes 1,2,..., one can 
define a permutation-ordered product as 

\d)p = cxp (ap^c'p^ - apjCPi^ cxp (ap^c'p^ - a*p^cp^'j ' ' ' |0>, (20) 
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which depends on P. For different permutations one obtains physically different states, so one has many different 
multi-mode "approximate complex fermion coherent states" (iV! variants for a finite number N of occupied modes). 
This fact is demonstrated in Fig.^ which shows how the complex degree of coherence (defined in AppendixQ depends 
sensitively on the choice of permutation P for two distinct multi-mode "approximate complex fermion coherent states" . 

(a) 



(b) 




FIG. 1: The square of modulus of the complex degree of coherence, |7(a;,y)p, for two "approximate complex fermion coherent 
states" according to Eq. 12011 that differ only by the permutation P. The modes are forty one-dimensional monochromatic 
waves with equally spaced wavenumbers ("comb") ordered according to their increasing wavenumber, and all the amplitudes 
Qffc are chosen to be equal to 0.166. The permutation in case (a) is identity, in case (b) a random permutation. The correlators 
clearly differ, which shows that the ordering of displacement operators is of large importance and different orderings result in 
physically different states. 

One could attempt to remove the permutation dependence of Eq. (|20|l by averaging it over all permutations of the 
modes: 

R = ^|d)p, (21) 
p 

up to a normalization factor. However, the resulting state has only zero- and single-particle components because all 
the multi-particle parts cancel due to the fermion anticommutation relations. As the state \d) contains no more than 
one fermion, it reduces to a single-mode "approximate complex fermion coherent state" . Thus we see that defining a 
multi-mode fermion coherent state using displacement operators requires an a priori mode ordering. Such an ordering, 
however, cannot be justified physically as none of the modes is preferred above the others. 
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Linear mode coupler: — Suppose that the most general single-mode fermion pure state /3o|0) + Pi\l) is mixed with 
the vacuum state via linear mode coupling. The resultant output state is 

/?o|0)i|0)2 + t/3i|l)i|0)2 + r/3i|0)i|l)2, (22) 

where the indices 1, 2 label the output modes and t and r is the complex transmissivity and reflectivity, respectively. 
Evidently expression (|22|l cannot be written as a product of a state of the first and second output ports, respectively, 
unless /3i = 0. For small /3i, which corresponds to the above definition of the "approximate complex fermion coherent 
state", one obtains almost a product state in H22(l . Hence, "approximate complex fermion coherent states" approx- 
imately satisfy the linear mode coupler condition. However, the condition is satisfied exactly only by the vacuum 
state, that is, in the limit of no particles. 

Summary: — Although the single-fermion "approximate complex coherent states" approximately satisfy the co- 
herence requirements given above (where, by single-fermion, we mean a superposition of no fermion and one fermion 
with the one- fermion contribution very small), we shall see that multi-mode "approximate complex fermion coherent 
states" are not so accommodating. In the next section we state several important propositions that hold for normally- 
ordered fermion correlators and do not have direct counterparts for boson fields. Based on these results, we will show 
that the correlator factorization condition cannot be satisfied for "approximate complex fermion coherent states" , 
even approximately, except for the case of single-particle fermion states. 

V. PROPERTIES OF NORMALLY-ORDERED FERMION CORRELATORS 

Fermion correlators are defined analogously to the boson counterparts and provide information on expected mea- 
surement outcomes by instruments at different spacetime coordinates. Definitions and examples of fermion correlators 
are presented in Appendix lO 

In the following we prove several important propositions that hold for normally-ordered fermion correlators of an 
arbitrary state (pure or mixed) and that are consequences of Pauli's exclusion principle. 

Proposition 2. For any fermion field state, the fermion correlator F'"^ (xi , . . . , a;„, y„, . . . , i/i) tends to zero smoothly 
whenever two points Xi, Xj or ?/,;, yj approach each other. 

Proof. Without lost of generality we can assume that i — l,j ~ 2. Using the mode decomposition of tj}{xi) — 
ipkixi)ak, we obtain 

r^"^ = E 'pUxiWlix2){aial4'Hx3) ■ ■ ■ i^HxnWvn) ■ ■ ■ ^{yi)) 

kl 

= I Yl[^kixi)Vi{x2) - iPi{xi)ipl{x2)]{alaj4!''{x3) ■ ■ ■ {xn)i'{yn) ' ' ' '^(yi)), (23) 

kl 

where the anticommutation relation aj^.aj^ = — a|a^ was employed. 

In this way, if xi — > X2, the antisymmetric function (p^.(.Ti)(pj*(x2) — ipi{xi)ipl.{x2) goes to zero and so does the 
whole correlator. The same argument can be used if some pair of points j/,;, yj goes together. The smoothness follows 
from the smoothness of the mode functions ipk{x). □ 

Proposition [3 shows that it is impossible to find two fermions at points too close to each other, which is a manifes- 
tation of Pauli's exclusion principle. 

Proposition 3. // the fermion correlator F^"-* (xi, . . . , 2/„, . . . , yi), for n > 1, factorizes for some state of the 
fermion field, then F("^(xi, . . . , a;„, y„, . . . , i/i) is identically equal to zero. 

Proof. The factorization of the correlator implies that 

f(")(xi, ... ,x„, ?;„,... ,?/i) f{xi)---f{xn)g{yi)---g{yn)- (24) 

However, if Xi = Xj {i ^ j), then it follows from Proposition |21 that F'-"-'(a;i, . . . , x„, yn, ■ ■ ■ , yi) = 0. Thus some of the 
functions f{x),g{x) must be equal zero; hence F*^"^ (xi, x„, ?;„,..., ?/i) vanishes. □ 

Proposition 4. Let |7(a:;,y)| — 1 and hence \T{x, y)\'^ — T{x,x)T{y,y) for some state j?/') and a pair of points x,y. 
Then F^") (a;, y, xa, . . . , x„, y„, . . . , yi) = for any n> 1 and any X3, . . . ,x„,yi, . . . ,y„. 
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FIG. 2: The second-order correlator [V'''^'' {x,y,y,x)\^ for a weak "approximate complex fermion coherent state" as in Fig.Q(a). 
The dip along the line a: = y is a common feature of this correlator regardless of the state and is a consequence of the fact that 
r*-^-* {x, y, y,x) ^ when x ^ y. 



Proof. The proof is based on argument similar to the one used in [lol.llll| for bosons in the case of the electromagnetic 
field. Define the operator C 

c = I i^{x) / ■■'^{y)- (25) 

Thus (^IC'I'CIV') = 1 - \l{x,y)\'^, which is equal to zero by assumption. It follows that C\ili) = 0. Hence 

This equation then yields 'il}{x)'ip{y)\'ij}) ~ because [■(/'(y)]^ = due to anticommutation relations. Hence, any 
normally-ordered correlator containing the operator product ?/'(x) ■(/;(?/) or ■(/;^(a;)'0t(j/) turns into zero. In particular 
r(")(x,y,...) = 0. □ 

One of the consequences of Proposition 0] is that whenever the field is perfectly first-order-coherent at the points 
x^y, it is impossible to find a fermion at the point x and another fermion at y. This can be understood as another 
(perhaps unusual) demonstration of the Pauli exclusion principle. 

Proposition 5. Let \'y{x,y)\ = 1 for some state \'ip) for all x,y. Then the state has support over only the vacuum 
and single-particle states. That is, the probability for two fermions in the field is identically zero. 

Proof. The proof follows directly from Proposition^ An alternative proof is based on an related result for optical 
field stating that if |7(a;,j/)| — 1 Vx, y, then only one mode of the field is occupied This holds also for fermions 
and in combination with Pauli's principle it yields that there is at the most one fermion in the field. □ 

Impossibility of factorizing fermion correlators: — The above propositions, and Proposition |3 in particular, show 
that the fermion correlators of order larger than one cannot factorize except for the vacuum and single-particle 
states. This expresses a deep difference between boson and fermion fields: whereas the requirement of the Glauber 
factorization leads to coherent states for bosons, for fermions it leads to states containing zero or one particle. 

Proposition |S1 imposes another strong condition on the fermion field. In particular, a field with more than one 
fermion cannot be perfectly first-order-coherent for all x,y. Indeed, if for some field, r<-^\x,y) = f{x)g{y) holds for 
all X, y, then 17(^^(0;, y)\ = 1. Consequently, according to Proposition|Sl it follows that the field does not contain more 
than one fermion. Hence, even the first-order correlators cannot factorize for most fields. Importantly, even first-order 
correlators cannot factorize for any fermion field with support over multi-particle (even two-particle) states. 
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First-order correlator for "approximate complex fermion coherent states": — To demonstrate the result from the 
previous section, we will now calculate the first-order correlator for the multi-mode "approximate complex fermion 
coherent state" (|2n|l with P the identity permutation: P = id. Using the decomposition aj ~ |Q!j|e"^^ we get: 

r^'Hx,y) =idH^^(x)7A(2/)|a)id 

= \ E ^*i^)^i(y) sin(2|a,|) sin(2|a, I) U,, e'(<^^-*') + ^ ^*ix)if,iy) sin' |a,|. (27) 

Here Uij is the product of the factors cos(2|Q;fc|) with k taking all values between i and j (with z, j themselves excluded). 
Hence, for i -I- 1 < j 

C/y =cos(2|a,+i|)cos(2|a,+2|)---cos(2|aj_i|), (28) 

and for j + 1 < i we use Eq. (|28|l with i,i interchanged. The factor Uij comes from the anticommutation relations 
between the field operators. 

If it were not for the factor Uij in Eq. (|27|l . then the correlator would approximately factorize if all \ai\ were 
small. However, the factor Uij causes the correlator to be far from factorization even for small \ai\. This can be 
seen in Fig. ^ if T^^\x,y) factorized, then the complex degree of coherence \"f'^^\x ^ y)\^ would be equal to unity. 
Clearly, this does not happen in Fig. ^ although the mean occupation number of each mode is very low: specifically, 
(sin|a|)^ = 0.027. Comparing this situation with bosons where the correlation function G^^\x,y) factorizes for a 
multi-mode coherent state, we see that there is a fundamental difference between fermions and bosons even in the 
limit of very low occupation numbers per mode. 

The above results show that higher-order correlators cannot be factorized nontrivially for any state of the fermion 
field and that first-order correlators can be factorized approximately only for states containing no more than one 
fermion. In the case that there is truly zero or one particle, we expect the features of fermions to be the same as 
those of bosons because at least two particles are required to see manifestations of their quantum statistics. However, 
higher-order correlators that describe two- or multi-particle characteristics of the field do reveal the fundamental 
difference between fermions and bosons regardless of how weak the field is because the correlator tells us about the 
correlation statistics regardless of how rare the multi-particle events are. This is further demonstrated in Appendix IDl 
on an example of boson and fermion chaotic states. 

Hence we see that there are no fermion states with complex parametrization that would, at least approximately, 
factorize fermion correlators, with the exception of single-particle states. Thus the definition of coherent states in 
terms of correlator factorization does not yield reasonable approximate multi-particle fermion coherent states, similarly 
as the previous definitions. 

VI. CONCLUSION 

Constructions of complex fermion coherent states, or approximations to such states, is problematic, as we have shown 
in this paper. We have seen that multi-mode "approximate coherent states" , whether obtained as displaced vacua or 
as near-eigenstates of the annihilation operator, exhibit inconsistencies due to inequivalences under permutation of 
modes, and also do not factorize fermion correlation functions. The only fermion states that approximately satisfy 
the coherent state definitions are the single-mode "approximate coherent states" that contain, however, no more than 
one particle. 

These problems arise ultimately from Pauli's exclusion principle and the corresponding mathematical framework 
of anticommuting operators. The ambiguity of mode ordering, the complications of linear mode coupling, and the 
unfactorizability of the correlation functions are due to the exclusion principle. On the other hand, for states that 
contain no more than one particle, the Pauli principle has no effect and "approximate coherent states" can be defined 
in the subspace of the Hilbert space spanned by zero- and single-particle states. In this domain of the Hilbert space, 
there is no essential difference between fermions and bosons as single particles cannot manifest their distinct quantum 
statistics. Hence single-mode "approximate complex fermion coherent states" approximately satisfy the coherence 
properties in the same way as "approximate boson coherent states" , which are obtained from weak boson coherent 
states by truncating the particle number beyond one. 

It therefore appears that the use of Grassmann numbers, which obviate the problem of anticommuting operators 
but introduce an unphysical parametrization, are needed to create self-consistent fermion coherent states that satisfy 
the expected coherence properties 0; complex parametrized fermion coherent states do not appear to be useful, even 
approximately. 
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APPENDIX A: MODE ORDERING FOR MIXTURES OF NUMBER STATES 

A mixture of fermion number states does not suffer from the permutation ambiguity described in Sec.^ Eqs. © 
and (|10|l . Indeed, the projector to a number state is invariant with respect to changes of the mode ordering due to 
the doubled number of the minus signs arising from the field operator anticommutation relations. For example, 

-P(iii2) = |lil2>(l2li| = aI4|0i02)(020i|a2ai (Al) 
A12I1) = |l2li)(lil2| = 4ai|0i02)(020i|aia2 (A2) 
= (A3) 

It follows that also an arbitrary mixture of number state projectors, e.g. a chaotic state, is immune with respect to 
the change in ordering the modes. This is a contrast to superpositions of number states, for which the ambiguity is 
present and changing the mode ordering can lead to a physically different state, as is discussed in Sec. IIVI 



APPENDIX B: PROOF OF EQUIVALENCE OF DEFINITIONS OF BOSON COHERENT STATE 

Here we prove that all the four definitions of boson coherent state mentioned in Sec. Illll are equivalent. We denote 
the properties as follows, 

(a) boson coherent state is an eigenstatc of the annihilation operator, 

(b) boson coherent state is a displaced vacuum state \a) = e"" ~" ° |0) = D{a) |0) 

(c) boson coherent state is a pure state that yields a product state when mixed with the vacuum state on a linear 
mode coupler (linear mode coupler). 

(d) boson coherent state is a state for which all normally-ordered correlators factorize. 

We will show now the implications (a) (c), (c) ^ (b), (b) ^ (a), and (a) (d). To show the implication (d) 
(a) is not so easy and for the proof we refer to the paper 

First we will show (a) (c). Let us denote the annihilation operators of the input and output ports of the linear 
mode coupler by ai, 0,2 and a'l, a'2, respectively. The relation between cii, 0,2 and a'l, a'2 is 

2 

with Aij a unitary matrix. Let an eigenstate \a) of ai incide on the first linear mode coupler input port, and vacuum 
incide on the second input port, which is an eigenstate of 0,2 with eigenvalue zero. It then follows from Eq. HBlfl that 
the linear mode coupler output state is an eigenstate of a'l and 03, which is a product state of the linear mode coupler 
outputs. 

Next we show (c) ^ (b): Let the boson coherent state incident on linear mode coupler be expressed in the Fock 
basis as follows: 

00 

|a)=5]c„(at)"|0) (B2) 

n—O 

The series X]^o'^"(^^)" formally defines a function of the creation operator that we will denote by f{a^)- The 
linear mode coupler input state is then f{a\)\0) as there is vacuum at the second input port. The relation between 
a\,a2 ^'^d a[',ay following from Eq. (|B1|1 is a] = X]j=i ^ji^'' Therefore the output linear mode coupler state is 
f{Aiia\' + Ai2a\')\0) . According to our assumption, this is a product state of output modes 1', 2', so it must hold 

f{Ana\' + Ai24') = g{a\')h{al') (B3) 
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with g, h some functions. As a{' and al' commute, we can treat /, g, h as ordinary functions of, say, arguments x and 
y instead of a\ and a\. Then, denoting F{z) = ln/(z). making logarithm of Eq. I|B3|) and performing the derivative 
d'^/dxdy, we arrive at the condition 

AnAi2F"{Anx + Ai2y) = 0. (B4) 

For both An and A12 different from zero (the opposite would correspond to a trivial linear mode coupler — either 
doing nothing or interchanging the modes 1 and 2), this is possible only if _F is a linear function, i.e., if F{z) = Az + 7. 
It then follows that /(fit) = e^"^+''. However, then fia\)\0) is a displaced vacuum state (up to a normalization 
constant) because the CampbcU-Hausdorff-Bakcr formula yields 

gAat+T,|Q^ ^ ^,|A|V2+7c^at-A*a|Q^, ^ e^^^^ D{X)\0) . (B5) 

To show (b) => (a) means to show that D{a)\0) is an eigenstate of the annihilation operator. Displacing the 
annihilation operator, we have 

aD{a)\0) = D{a)&{a)aD{a)\0) ^ D{a){a + al)|0) = aL»(a)|0), (B6) 

which we wanted to prove. 

Next we show (a) =^ (d): The field operators and -0^(2;) are linear combinations of single-mode annihilation 
operators and creation operators d|, respectively. Any (possibly multi-mode) coherent state is therefore a right 
eigenstate of V'(x) and left eigenstate oii>'^{x). From this the factorization H16|l follows immediately. 



APPENDIX C: NORMALLY-ORDERED CORRELATORS 



Correlators are important tools for describing coherence of both boson and fermion fields. They are defined as 
expectation values of products of creation and annihilation operators at different space-time points. 
The fermion correlators of the n"^ order are defined by analogy to their boson counterparts 

r^")(xi, ... ,a;„, ?/„,... ,yi) = (V''^(a;i) • • • Vi^(a;„)'0(yn) • ■■^{yi)) (CI) 

(we use the Greek letter F to distinguish from the boson correlator G). Of particular importance is the correlator 
with repeated arguments, for which we introduce the shorter expression F'"-'(xi, . . . ,a;„) = F'"-'(a;i, . . . ,x„,a;„, . . . ,a:i), 
which corresponds to the ri-fermion detection probability, and also the one-fcrmion cross-correlator F(x, y) = F'^' (x, y). 

The normalized first-order cross-correlator '~f{x,y) = T{x, y) / ^jT{x, x)T{y, y) is called complex degree of coherence 
and its magnitude expresses the visibility of interference fringes in a double-slit experiment in which the two slits are 
placed at the points x and y, respectively If 2/)l = 1 for 2;, y, we say that the field has perfect first-order 
coherence. 



APPENDIX D: CORRELATORS FOR BOSON AND FERMION CHAOTIC STATES 



We will show now that for chaotic states, the fermion or boson nature of the particles does not have an effect on 
the first-order correlator while it has a strong effect on the higher-order correlators. 

For a chaotic state |l5l | , each mode of the chaotic field has maximum entropy for a given mean number of particles 
in this mode, and the density operator of the field is a direct product of single-mode density operators. This second 
property simply implies that, in the chaotic state, the individual modes of the field are totally uncorrelated. 

The chaotic state density operator of a single mode of bosons or fermions with mean number of particles M can be 
expressed as 



with 



'^^iTm' Tr(e-«"'") = 1 + M (D2) 
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for bosons and 

e"^ = T^n^ ^ (e-«°'°) = —1— (D3) 

for fermions (unlike previous sections, here we denote the fermion field operators by a, a) rather than c, for compact 
notation). Thus we have the density operators for bosons and fermions in the following forms: 



PB = 



n=0 ^ ' n=0 ^ ^ 



Consider now a multi-mode chaotic field of bosons or fermions. Let there be N occupied modes labeled by 1, . . . , A'' 
and let the mean number of particles in the i*'^ mode be Mi . For brevity, we denote the boson and fermion correlators 
by Gb and Gp, respectively, in this section. The first order cross-correlator Ggp(a;, y) can readily be calculated and 
has the same form for bosons and fermions: 

N 

1=1 

Eq. IjDSp shows that the first order coherence (which is related to the visibility of interference fringes) for chaotic 
states is not at all influenced by the fermion or boson nature of the particles, provided the mean occupation numbers 
are the same. This can be expected as the first order correlator is not related to multi-particle correlations so it should 
not be affected by the exchange interaction. 

The n*'^-ordcr correlator can also be calculated for the chaotic state. The following formula holds exactly for 
fermions and for bosons it holds with good accuracy if Mi ^ 1,« = 1,...,A^: 



G^"^(a;i, . . .,Xn,yn, ■ • ■ ,yi) = X! ^b?f^''^b,f(^1' yp<i))^B,F(^2, yp(2)) • ■ ■G'^^p{xn,yp(,^)) (D6) 

p 



Here the sum runs over all permutations P of the indices 1, 2, . . . , rt, the parity of the permutation P is denoted by 
par(P), and xb = 1, Xf = ^ 1 is the boson and fermion sign factor, respectively. Hence the higher-order correlators 
are no longer identical for boson and fermion chaotic states, which is manifested as bunching (antibunching) for boson 
(fermion) chaotic fields. If Xi Xj,i ^ j, Gf goes to zero according to Eq. (|D6|I due to the factor (^—l'jP^'-'(P\ which 
is in agreement with the general property of fermion correlators (Proposition [5J). 
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